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1.0  SUMMARY 

We  have  examined  some  of  the  methods  available  for  use  in  the  labora¬ 
tory  to  measure  the  inelastic  attenuation  of  seismic  waves  in  rocks.  We  con¬ 
clude  that  the  forced  swept  resonance  technique  is  best  suited  for  measuring  the 
Q  of  linear  anelastic  materials,  and  for  defining  the  transition  amplitude  from 
linear  anelastic  behavior  to  nonlinear  behavior.  For  evaluating  losses  in 
materials  at  high  amplitudes,  which  exceed  the  elastic  limit,  the  best  technique 
requires  simultaneous  measurements  of  the  time  functions  of  stress  and  strain 
under  conditions  that  simulate  seismic  loading  as  closely  as  possible. 

We  have  also  studied  the  transition  region  from  linearity  to  nonlinear¬ 
ity  in  a  test  specimen  of  Westerly  granite  at  elevated  effective  pressures  using 
the  mechanical  resonance  approach.  Nonlinear  effects,  both  in  flexure  and  in 
shear,  are  observed  when  the  strain  exceeds  10~6,  increasing  slightly  with  in¬ 
creasing  effective  pressure.  These  transition  amplitudes  probably  represent 
lower  limits  on  the  amplitude  of  transition  from  linear  to  nonlinear  behavior 
for  a  compressi onal  pulse  propagating  through  the  near-fields  of  an  explosion. 

The  details  of  stress-strain  hysteresis  loops  have  been  examined  for  a 
number  of  different  rock  specimens  using  alternating  compressive  and  tensile 
stresses.  Our  experimental  results  support  the  following  general  observations: 
(1)  rocks  have  a  larger  effective  modulus  in  compression  than  in  tension; 

(?)  most  of  the  energy  lost  during  a  full  cycle  of  loading  occurs  as  a  result  of 
strain  in  extension,  since  the  hysteresis  loops  are  smaller  in  compression  than 
in  tension;  (3)  stress-strain  hysteresis  loop  traces  appear  to  be  independent  of 
frequency.  These  observations  indicate  a  loss  mechanism  associated  with  inter¬ 
granular  friction.  In  compression,  intergranular  sliding  is  apparently  restric¬ 
ted  by  the  impingement  of  opposing  crack  faces,  while  in  tension  this  does  not 
occur.  Thus,  laboratory  measurements  of  attenuation  that  use  alternating  com¬ 
pressive  and  tensile  loads  may  overestimate  the  attenuation  of  high  amplitude 
compressive  waves  with  the  same  maximum  strain. 
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2.0  INTRODUCTION 

The  mechanics  of  wave  propagation  through  the  near  field  of  an  explo¬ 
sive  seismic  source  are  not  completely  understood,  particularly  in  the  peak 
strain  amplitude  range  between  10"^  and  10“^.  The  most  complete  data  sets  that 
describe  free-field  wave  propagation  as  a  function  of  scaled  distance  (distance/ 
yield^/3)  in  this  strain  amplitude  range  have  been  obtained  for  explosions  in 
polycrystalline  salt.  The  data  base  includes  measurements  from  (a)  the  5.3  kT 
nuclear  event,  SALMON  (Perret,  1967;  Trulio,  1978);  (b)  the  medium  to  large 
scale  chemical  explosions  of  the  COWBOY  experiments,  which  took  place  in  a 
natural  salt  done  (Trulio,  1978,  1931);  and  (c)  a  series  of  small  scale  chemical 
explosions  in  pressed  polycrystall ine  salt  in  the  laboratory  (Larson,  1982). 

When  all  available  data  are  combined,  it  appears  that  wave  propagation  satisfies 
cube  root  of  yield  scaling.  That  is,  the  decay  of  peak  particle  velocity  and 
displacement  can  be  defined  as  a  function  of  distance/yield1/^  (Trulio,  1978, 
1931;  Larson,  1982).  This  appears  to  be  the  case,  even  though  available  data 
cover  approximately  10  decades  in  yield,  4  decades  in  peak  particle  velocity, 
and  4  decades  in  frequency.  However,  the  data  fit  a  line  that  decays  much 
faster  than  r"*,  which  should  be  the  case  if  the  material  behaves  in  a  perfectly 
elastic  manner.  Thus,  the  behavior  of  salt  in  the  range  of  available  free  field 
measurements,  to  strains  as  low  as  approximately  7  *  10“^,  cannot  be  regarded  as 
perfectly  elastic  (Trulio,  1978,  1981).  Because  of  the  uncertainties  regarding 
the  amount  of  erergy  that  can  be  dissipated  as  a  seismic  pulse  propagates  at 
intermediate  strains,  Bache,  et  al  (1981)  have  questioned  the  usefulness  of 
reduced  displacement  potential  (RDP)  calculations  based  on  close  range  data  for 
the  purpose  of  defining  a  seismic  source  function. 

The  exact  nature  of  attenuation  in  the  near  field  of  an  explosive 
seismic  source  remains  a  controversial  subject.  In  fact,  the  evidence  regarding 
the  issue  of  linearity  vs  nonlinearity  appears  to  be  a  bit  ambiguous.  Larson 
(1982)  measured  a  Q  of  12.5  near  10"^  strain,  and  a  0  of  24.9  near  6  *  10-4 
strain,  in  the  laboratory  using  small  scale  chemical  explosions  in  pressed  salt. 
The  fact  that  0  is  amplitude  dependent  is  evidence  for  nonlinear  response.  How¬ 
ever,  Larson  also  demonstrated  an  approximately  linear  superposition  of  wave- 
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forms  at  strains  higher  than  10"4,  which  suggests  near-linearity.  The  experi¬ 
mental  results  reported  in  Tittmann  (1983)  using  resonating  bars  of  natural 
halite  indicate  that  the  linear  anelastic  Q  of  halite  is  quite  high,  near  500, 
with  a  transition  to  an  amplitude  dependent  nonlinear  Q  at  strain  amplitudes 
greater  than  approximately  2  *  10"®.  Burdick,  et  al.  (1984a)  have  argued  that 
it  is  possible  to  model  a  seismic  source  function  for  the  Amchitka  tests, 
detonated  in  volcanic  material,  assuming  linear  behavior  in  the  near  field  just 
outside  the  spall  zone  (~  700-1200  m/kT^/^).  Furthermore,  Burdick,  et  al 
(1984b)  contend  that  the  same  model  can  be  used  to  predict  the  first  vertical 
pulse  arrival  and  rise  times  even  within  the  spall  zone.  They  used  tne  concept 
of  a  congressional  elastic  radius,  which  in  fact  may  be  considerably  smaller 
than  a  tensional  elastic  radius,  and  which  must  extend  at  least  as  far  as  the 
outer  limits  of  the  spall  zone.  Minster  and  Day  (1985)  recently  re-examined  the 
COWBOY  data  set,  and  concluded  that  it  is  possible  to  explain  simultaneously  the 
radial  decay  of  peak  displacement  and  peak  velocity  by  either  (a)  a  linear 
anelastic  model  with  low  Q  (~  20)  or  (b)  a  nonlinear  model  with  amplitude- 
dependent  0.  MacCarter  and  Wortman  (1985)  conclude  that  the  free-field  motion 
measurements  from  the  SALMON  event  are  consistent  with  an  ampl l tude- i ndependent 
0  of  about  10.  In  any  case,  these  amplitude-independent  Q  values  are  too  low  to 
represent  the  linear  anelastic  Q  of  salt  according  to  the  results  of  Tittmann 
(1983). 

In  a  more  recent  study  of  explosive-induced  wave  propagation  in  Nugget 
sandstone,  Larson  (1984)  reported  Q  measurements  between  8  and  20  at  strains 
between  2  *  lO"4  and  2  *  10"^.  Although  low  strain  (<  10'®)  attenuation  meas¬ 
urements  are  not  available  for  Nugget  sandstone,  Tittmann  et  al  (1981)  and 
Winkler  et  al  (1979)  report  relatively  high  Q  values  (>  100)  for  dry  Berea  and 
Boise  sandstones  under  linear  conditions. 

While  generally  it  is  acknowledged  that  Q  is  defined  rigorously  only 
when  the  material  through  which  a  wave  propagates  behaves  linearly,  it  is  common 
(McCarter  and  Wortman,  1985;  Minster  and  Day,  1985)  to  assume  for  theoretical 
purposes  that  "nonlinear  Q"  can  be  defined  using  the  equations  of  Mavko  (1979) 
and  Stewart,  et  al  (1983): 
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Q  *  (w,e)  =  Qg'^(w)  +  ae  (1) 

where  Q"*  is  nonlinear  attenuation,  Qq"*  is  linear  anelastic  attenuation,  u  is 
frequency,  a  is  a  constant,  and  e  is  strain.  This  equation  considers  the  com¬ 
bined  effects  of  Coulomb-type  friction  from  many  intergranular  contacts  on  the 
energy  dissipated  during  one  full  elastic  wave  cycle.  In  their  derivations,  Q”* 
was  defined  as 

O"1  =  (1/2  u) (  AW/W)  (2) 

where  AW  is  the  energy  dissipated  in  one  full  cycle,  and  W  is  the  peak  strain 
energy  stored  per  cycle. 

One  technique  used  in  the  laboratory  to  estimate  nonlinear  Q  involves 
the  forced  resonance  of  bars.  The  resonating  bar  approach  offers  a  particularly 
sensitive  tool  for  detecting  the  onset  of  nonlinearity,  even  at  very  low  strains. 
The  data  reported  by  Tittmann  (1983)  indicate  that  the  intrinsic  linear  attenu¬ 
ation  of  dome  salt  is  quite  low  (Q  >  500),  and  that  the  onset  for  nonlinear  be¬ 
havior  is  near  2  *  10-6  strain,  corresponding  to  approximately  1  bar  of  stress. 
Similar  behavior  has  been  reported  for  other  rock  types,  including  granite, 
limestones,  sandstones,  and  miscellaneous  igneous  rocks  (cf  Mavko,  1979;  Stewart 
et  al  ,  1983).  These  measurements  suggest  that  nonlinear  response  should  persist 
to  large  scaled  distances  from  explosions,  on  the  order  of  10^  m/kt*'  .  Fur¬ 
thermore,  the  low  Q  values  calculated  by  Larson  (1982)  [laboratory  measure¬ 
ments],  Trulio  (1979,  1981)  [COWBOY],  Minster  and  Day  (1985)  [COWBOY]  and 
McCarter  and  Wortman  (1985)  [COWBOY,  SALMON]  from  free-field  measurements  in 
salt  are  much  lower.  This  constitutes  additional  evidence  of  nonlinearity  in 
available  free-field  measurements. 

The  experimental  results  obtained  in  the  previous  program  (Tittmann, 
1983)  indicate  that  natural  dome  salt  is  rather  ideally  behaved.  Linear 
anelastic  behavior  is  indicated  when  (1)  the  attenuation  and  resonant  frequency 
are  independent  of  vibration  amplitude  and  (2)  a  proportionality  exists  between 
driving  voltage  and  vibration  amplitude.  At  elevated  effective  pressures  linear 
anelastic  behavior  is  observed  at  strain  amplitudes  below  about  2  *  10"®,  corre- 
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sponding  to  stresses  of  about  1  bar.  The  linear  anelastic  Q  for  a  dry  natural 
dome  salt  specimen  was  quite  high  (~  500  in  extension  and  ~  1000  in  torsion), 
apparently  independent  of  pressure  up  to  6.8  x  10^  Pa,  and  very  slightly  de¬ 
pendent  on  frequency  between  ~  80  Hz  and  ~  480  Hz.  Nonlinear  behavior,  as  re¬ 
flected  in  an  amplitude-dependent  0,  was  observed  at  strain  amplitudes  above 
2  x  lCT®,  at  all  effective  pressures  to  6.8  x  io^  Pa  (corresponding  to  burial 
depths  of  ~  2.4  km).  The  behavior  of  pressed  salt  contrasted  with  that  of 
natural  dome  salt.  Weakly  nonlinear  behavior  persisted  to  very  low  amplitudes 
(below  10"®),  corresponding  to  stresses  of  approximately  0.01  bars.  Strongly 
nonlinear  behavior  is  observed  at  strain  amplitudes  higher  than  10"®.  Under 
ambient  pressure  conditions  attenuation  in  pressed  salt  was  sensitive  to  the 
amount  of  moisture  adsorbed  either  on  the  surface  or  within  the  interior  of  the 
specimen . 

The  previous  study  (Tittmann,  1983)  raised  a  number  of  questions  re¬ 
garding  the  experimental  methods  used  to  measure  nonlinear  0,  and  to  parameter¬ 
ize  the  nonlinear  inelastic  properties  of  rocks.  Solutions  to  some  of  these 
issues  were  explored  in  this  program  and  are  contained  in  this  report. 
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3.0  Q  AND  MODULUS  MEASUREMENT  TECHNIQUES 

3.1  Resonance  Methods 


3.1.1  Linear  Regime 

Several  of  the  issues  addressed  in  tnis  report  relate  to  using  resonat¬ 
ing  bar  measurements  for  estimating  P-  and  S-wave  velocity.  The  resonance  meas¬ 
urements  are  normally  carried  out  with  relatively  thin,  inertially  loaded  speci¬ 
mens  and  involve  0  and  modulus  measurements  corresponding  to  at  least  two  modes 
of  vibration:  flexure  or  extension,  and  torsion.  Calculations  of  seismic  body 
wave  properties  from  laboratory  measurements  are  reasonably  good  provided  that 
(1)  the  resonator  is  properly  designed  and  constructed,  and  (2)  the  response  of 
the  specimen  to  loading  is  macroscopical ly  homogeneous. 

The  apparatus  which  was  used  in  the  previous  program  (Tittmann,  1983) 
has  also  been  used  for  parts  of  this  study.  It  is  illustrated  schematically  in 
Fig.  3.1.  This  apparatus  is  used  to  generate  vibrations  and  to  measure  their 
amplitude  in  both  flexure  and  torsion.  Its  design  has  evolved  from  earlier 
models  to  include  the  following  features  which  are  essential  for  obtaining  reli¬ 
able  measurements  of  modulus  and  attenuation:  (1)  a  very  stiff  cl  amp-arm-rotor 
assembly  designed  to  minimize  bending  moments  in  any  part  of  the  apparatus  other 
than  the  test  specimen,  (2)  rigid  mechanical  coupling  between  the  test  sample 
and  the  rest  of  the  resonator  assembly,  and  (3)  minimal  internal  friction  within 
the  apparatus  itself  during  sample  resonance.  The  resonator  has  been  designed 
to  minimize  instrumental  corrections.  However,  careful  measurements  have  been 
performed  using  low  dissipation  standards  with  known  modulus  in  order  to 
quantify  instrumental /system  factors,  such  as  aerodynamic  drag,  which  must  be 
incorporated  into  modulus  and  attenuation  calculations. 

3.1.2  Nonlinear  Regime 

Although  the  resonant  bar-type  measurements  are  well  suited  for  measur¬ 
ing  attenuation  and  modulus  in  the  linear  anelastic  regime,  and  for  defining  the 
transition  from  linear  to  nonlinear  behavior,  several  other  issues  need  to  be 
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Fig.  3.1  Schematic  illustration  of  the  resonant  bar  apparatus  used  for  studies 
of  dynamic  modulus  and  attenuation  as  a  function  of  strain  amplitude. 


considered  before  one  can  apply  the  high  amplitude  (nonlinear)  Q  measurements  to 
studies  of  high  amplitude  seismic  pulse  propagation  in  the  near  field. 

1.  Forced  Resonant  Vibrations  -  Q  can  be  calculated  using  the  formula 


where  fp  is  the  resonant  frequency,  measured  at  peak  amplitude,  and  Af  is  the 
bandwidth  of  the  peak,  measured  at  half-power  amplitude.  It  is  well  established 
(cf.  Nowick  and  Berry,  1972)  that  this  is  a  sound  technique  for  measuring  0  in  a 
material  with  linear  response.  However,  this  may  not  be  a  particularly  satis¬ 
factory  technique  for  measuring  energy  dissipation  in  a  material  with  properties 
that  depend  on  vibration  amplitude.  In  most  rocks,  when  vibrations  exceed  the 
elastic  limit  the  modulus  decreases  with  increasing  vibration  amplitude  (Mavko, 
1979;  Tittmann,  1983).  This  is  probably  due  to  frictional  sliding  along  micro- 
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cracks,  a  nonlinear  and  nonrecoverable  process.  The  shape  of  the  resonance  peak 
at  low  amplitude  is  nearly  symmetrical.  At  high  amplitudes  the  peak  becomes 
severely  skewed  toward  the  high  frequency  side.  This  effect,  illustrated  in 
Fig.  3.2  for  natural  dome  salt  at  0.68  MPa  effective  pressure,  is  a  result  of 
the  decrease  in  modulus  with  increas¬ 
ing  vibration  amplitude,  and  has  been 
observed  at  high  amplitudes  in  all  but 
the  most  tightly  consolidated  rock 
types.  This  result  casts  suspicion  on 
the  usefulness  of  the  swept  forced 
resonance  technique  for  obtaining  re¬ 
liable  measurements  of  Q  and  modulus 
at  high  nonlinear  amplitudes. 
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Fig.  3.2  I 

Resonance  curve  shapes  at  various 
strain  amplitudes.  Low  amplitude 
vibrations  within  the  linear  regime 
are  generally  nearly  symmetric,  as 
shown  in  the  lower  frame.  Skewness  is 
apparent  at  higher  nonlinear  ampli¬ 
tudes  (center  and  top  frames). 
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2.  Homogeneity  of  Strain  -  Another  problem  with  measuring  nonlinear  0 
involves  the  variation  in  strain  amplitude  throughout  the  speci¬ 
men.  In  the  previous  study  (Tittmann,  1983)  we  reported  measure¬ 
ments  of  Q  as  a  function  of  strain  amplitude,  which  was  calculated 
at  a  position  of  maximum  strain  in  the  specimen.  In  these  meas¬ 
urements  first-order  strains  varied  from  zero  to  the  value  actu¬ 
ally  reported  for  both  torsional  and  flexural  vibrations,  depend- 
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ing  on  location  within  the  specimen.  Maximum  strain  amplitude  is 
defined  for  torsion  and  flexure  in  Figs.  3.3  and  3.4,  respectively. 
Minster  and  Day  (1985)  discuss  this  issue  in  more  detail  and  esti¬ 
mate  the  magnitude  of  its  effect  on  the  attenuation-amplitude 
curves . 


•ci*  2i*a< 

TORSION  STRAIN  AMPLITUDE 


tea. 

FLEXURE  STRAIN  AMPLITUDE 


Fig.  3.3  Definition  of  strain 
amplitude  in  torsion. 


Fig.  3.4  Definition  of  strain  am¬ 
plitude  in  flexure. 


3.  Tension  vs  Compression  -  During  both  extensional  and  flexural 
resonances  each  increment  of  volume  within  the  specimen  is  sub¬ 
jected  to  tensile  stresses  half  of  the  time  and  to  compressive 
stresses  the  other  half.  A  tacit  assumption  is  that  the  material 
behaves  the  same  in  compression  as  it  does  in  tension.  The 
apparent  drop  in  bar  modulus  with  increasing  vibration  amplitude 
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has  been  noted  previously  (Tittmann,  1983),  and  is  probably  also 
related  to  the  distortion  of  the  resonance  peak  shape  with  in¬ 
creasing  strain.  It  has  not  been  proven,  however,  whether  either 
the  frictional  losses  or  the  effective  modulus  of  a  rock  are  the 
same  in  tension  as  they  are  in  compression.  Resonating  bar-type 
measurements  are  not  capable  of  resolving  the  details  of  the 
rheological  response  of  materials  to  nonlinear  loading. 

3.2  Quasi-Static  Stress-Strain  Analysis 

On  the  basis  of  the  above  discussion  it  is  evident  that  the  laboratory 
work  most  beneficial  to  resolving  the  issue  of  linearity  vs  nonlinearity  would 
involve  studies  of  the  mechanical  response  of  various  test  site  materials,  in¬ 
cluding  salt,  to  realistic  seismic  pulse  loading  histories.  Such  studies  would 
necessarily  require  the  application  of  true  umaxial  strain  while  test  specimens 
are  subjected  simultaneously  to  a  hydrostatic  bias  stress.  This  type  of  study 
is  very  difficult  technically,  and  results  of  such  studies  do  not  exist,  at 
least  not  in  the  open  literature.  Several  papers  have  examined  the  details  of 
stress-strain  hysteresis  loops  under  nonlinear,  uniaxial  stress  loading 
conditions  (distinguished  from  uniaxial  strain)  (Gordon  and  Davis,  1968; 
McKavanagh  and  Stacey,  1974),  but  in  each  of  these  cases  only  compressi onal 
measurements  were  performed. 

In  this  program  we  have  examined  the  effects  of  combined  high  amplitude 
tensile  and  compressive  uniaxial  loading  stress  on  four  different  rock  types: 
Westerly  granite,  Boise  sandstone,  Berea  sandstone,  and  Indiana  limestone,  under 
ambient  pressure  conditions.  This  represents  the  first  stage  of  a  new  effort 
with  the  ultimate  goal  to  examine  the  mechanical  response  of  test  site  materi¬ 
als,  including  dome  salt,  to  high  amplitude  uniaxial  strain  loading  at  elevated 
hydrostatic  pressures.  The  details  of  the  stress-strain  hysteresis  loops  are 
studied  in  order  to  determine  whether  loss  is  related  to  intergranul ar  friction 
or  to  linear  anelastic  relaxation.  The  response  of  the  rock  to  compressive 
loading  is  compared  with  that  of  tensile  loading. 
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4.0  MECHANICAL  RESONANCE  STUDIES 

We  have  examined  the  response  of  Westerly  granite  to  sinesoidal  loading 
using  flexural  and  torsional  resonances.  The  measurements  clearly  define  the 
amplitude  of  transition  from  linearity  to  nonlinearity  at  elevated  effective 
pressures.  The  results  are  compared  with  previous  measurements  of  dome  salt 
under  similar  conditions. 

4.1  Experimental  Techniques 

The  techniques  used  for  resonance  measurements  are  essentially  the  same 
as  t  cse  described  previously  (Tittmann,  1983).  Briefly,  a  cylindrical  specimen 
of  Westerly  granite,  13  cm  long  and  0.71  cm  radius,  was  clamped  into  an  appa- 
tus  used  to  force  flexural  and  torsional  vibrations.  In  this  particular  study 
the  resonant  frequency  was  between  500  and  550  Hz  for  all  measurements.  The  Q 
of  the  resonance  was  calculated  from  the  ratio  of  the  resonant  frequency  to  the 
bandwidth  of  the  resonance  curve.  The  voltage  applied  to  the  electromagnetic 
transducer  used  to  drive  the  resonance  was  varied  in  order  to  affect  the  ampli¬ 
tude  of  the  vibration.  Thus,  attenuation  and  RMS  driving  voltage  were  measured 
as  a  function  of  vibration  amplitude  in  the  sample.  Provided  that  the  sample 
behaves  linearly,  the  attenuation  should  be  independent  of  vibration  amplitude, 
and  the  RMS  driving  voltage  should  be  proportional  to  vibration  amplitude  (c.f, 
Tittmann,  1983).  Measurements  were  made  in  both  wet  and  dry  Westerly  granite  at 
effective  pressures  from  0  to  34  MPa. 

4.2  Experimental  Results 

The  results  of  this  study  are  shown  in  Figs.  4.1  through  4.20.  They 
are  qualitatively  consistent  with  the  results  of  previous  studies  of  natural 
dome  salt  (Tittmann,  1983)  in  that  linear  behavior  is  observed  at  all  low  ampli¬ 
tudes,  with  a  transition  to  nonlinear  behavior  at  high  strain.  The  transition 
amplitude  is  near  10"^  in  all  cases,  including  both  torsion  and  flexure,  regard¬ 
less  of  water  content,  increasing  very  slightly  with  increasing  effective 
pressure. 
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Fig.  4.1 

Attenuation  and  log  RMS  driving  trans¬ 
ducer  voltage  as  a  function  of  maximum 
strain  amplitude  for  torsional  vibra¬ 
tions  in  dry  Westerly  granite  at  0  MPa 
effective  pressure. 
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Fig.  4.2 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  func¬ 
tion  of  maximum  strain  amplitude 
for  torsional  vibrations  in  dry 
Westerly  granite  at  1.7  MPa 
effective  pressure. 
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Attenuation  and  log  RMS  driving  1 

transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
torsional  vibrations  in  dry 
Westerly  granite  at  6.8  MPa 
effective  pressure. 
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Fig.  4.5 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function  of 
maximum  strain  amplitude  for  tor¬ 
sional  vibrations  in  dry  Westerly 
granite  at  34  MPa  effective 
pressure. 
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Attenuation  and  log  RMS  driving  5 
transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
flexural  vibrations  in  dry 
Westerly  granite  at  0  MPa 
effective  pressure. 
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Fig.  4.7 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function  of 
maximum  strain  amplitude  for  flex¬ 
ural  vibrations  in  dry  Westerly 
granite  at  1.7  MPa  effective 
pressure . 
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Atter  .tion  and  log  RMS  driving 
transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
flexural  vibrations  in  dry 
Westerly  granite  at  3.4  MPa 
effective  pressure. 
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Fig.  4.9 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function  of 
maximum  strain  amplitude  for 
flexural  vibrations  in  dry  Westerly 
granite  at  6.8  MPa  effective 
pressure . 
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Fig.  4.10 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
flexural  vibrations  in  wet 
Westerly  granite  at  34  MPa 
effective  pressure. 
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Fig.  4.12 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
torsional  vibrations  in  wet 
Westerly  granite  at  1.7  MPa 
effective  pressure. 
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Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function  of 
maximum  strain  amplitude  for  tor¬ 
sional  vibrations  in  wet  Westerly 
granite  at  0  MPa  effective  pressure. 
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Fig.  4.13 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function  of 
maximum  strain  amplitude  for  tor¬ 
sional  vibrations  in  wet  Westerly 
granite  at  3.4  MPa  effective 
pressure . 
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Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
torsional  vibrations  in  wet 
Westerly  granite  at  6.8  MPa 
effective  pressure. 
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Fig.  4.15 

j  Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function  of 
maximum  strain  amplitude  for  tor¬ 
sional  vibrations  in  wet  Westerly 
^  granite  at  34  MPa  effective 
pressure. 
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Fig.  4.16 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function 
f  maximum  strain  amplitude  for 
exural  vibrations  in  wet 
esterly  granite  at  0  MPa 
ffective  pressure. 
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Fig.  4.18 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
flexural  vibrations  in  wet 
Westerly  granite  at  3.4  MPa 
effective  pressure. 
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Fig.  4.19 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function  of 
maximum  strain  amplitude  for 
flexural  vibrations  in  wet  Westerly 
granite  at  6.8  MPa  effective 
pressure. 
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Fig.  4.20 

Attenuation  and  log  RMS  driving 
transducer  voltage  as  a  function 
of  maximum  strain  amplitude  for 
flexural  vibrations  in  wet 
Westerly  granite  at  34  MPa 
effective  pressure. 
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5.0  QUASI-STATIC  STRESS-STRAIN  ANALYSIS 


5.1 


Experimental  Techniques 


(») 

LOAD  CELL 
(STRESS! 

HYDRAULIC 

GRIPS 

HYDRAULIC 

RAM 


MTS  LOAD  FRAME 


SAMPLE 


Figure  5.1  includes  schematic  il 1 ustrations  of  the  mechanical  parts  of 
the  system  used  for  measurement.  Briefly,  cylindrical  rock  specimens  2.54  cm  in 
diameter  and  6.35  cm  long  were  cut  from 
larger  blocks  of  material.  A  resistance- 
type  strain  gauge  was  bonded  to  the  sur¬ 
face  of  each  specimen  to  respond  to 
strains  in  the  longitudinal  direction. 

The  rock  was  then  bonded  adhesively  to 
stainless  steel  mounting  fixtures  using  a 
semi  plastic  polymer  resin  (Crystal  Bond 
^509,  Aremco  Products)  that  is  quite 
brittle  and  stiff  at  room  temperature,  but 
that  softens  readily  at  approximately 
70°C.  By  keeping  the  specimen  warm  while 
mounting  it  in  the  apparatus,  residual 
stresses  could  be  removed.  Measurements 
were  made  at  room  temperature  and  under 
ambient  pressure  conditions. 


EXTENSION  BARS 


'ADHESIVE  BONDING 


Fig.  5.1  Schematic  illustration  of  the  apparatus 
used  for  the  measurement  of  stress- 
strain  response  curves,  (a)  MTS  electro- 
fa)  MTS  electro-hydraulic  load  frame  and 
(b)  detail  of  sample  and  end  fixtures. 


Experiments  were  run  using  an  MTS  electro-hydraulic,  closed-loop  load 
frame  equipped  with  the  hydraulic  hardware  and  controls  necessary  to  apply  al¬ 
ternating  tensile  and  compressive  loads.  The  instrumentation  used  for  the  meas¬ 
urements  is  illustrated  in  Fig.  5.2.  The  frequency  of  loading  was  usually  1  Hz, 
except  for  two  series  of  measurements  on  Westerly  granite  and  Berea  sandstone  at 
0.1  Hz.  A  series  of  eight-cycle  bursts  of  load-controlled  constant  strain  rate 
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Fig.  5.2  Schematic  illustration  of  instrumentation  used  to 
acquire  and  process  stress-strai n  response  curves. 

triangle  wave  loading  were  applied  to  each  specimen.  Measurements  of  load  and 
strain  were  collected  simultaneously,  digitized,  and  stored  in  the  computer. 

Representative  time  functions  of  stress  and  strain  are  illustrated  in 
Fig.  5.3.  Measurements  were  obtained  on  each  sample  using  bursts  of  increasing 
maximum  stress,  starting  with  +3.4  x  1U^  Pa  and  working  up  to  a  maximum  of  +3.4 
x  107  Pa,  unless  precluded  by  the  breakage  of  the  specimen.  As  many  as  5  bursts 
of  8  waves  each  were  applied  at  low  amplitudes  in  order  to  enable  satisfactory 
signal  averaging  when  the  signal  to  noise  ratio  was  low.  At  higher  amplitudes 
the  si gnal -to-noi se  ratio  was  good,  and  only  one  eight-cycle  burst  was  run  in 
order  to  avoid  excessive  damage  to  the  specimen. 

5.2  Experimental  Results 

To  test  (a)  the  linearity  of  the  experimental  apparatus  and  (b)  the 
stiffness  of  the  adhesive  bonds,  the  first  phase  of  the  experimental  study  in¬ 
volved  measurements  of  stress  and  strain  on  an  aluminum  bar  which  has  a  rela¬ 
tively  ideal  linear  elastic  response.  A  typical  curve  of  stress  vs  strain  for 
an  aluminum  bar  is  shown  in  Fig.  5.4.  Extensional  stresses  and  strains  are 
positive.  Only  a  small  amount  of  hysteresis  is  observed  in  this  set  of  measure¬ 
ments,  and  it  is  apparent  that  the  relationship  between  stress  and  strain  is 
nearly  linear.  A  small  amount  of  hysteresis  is  observed  and  may  be  attributed 
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to  plastic  flow  in  the  resin  used  to  bond  the  aluminum  bar  to  the  loading  fix¬ 
tures.  This  is  a  relatively  insignificant  effect  compared  with  the  large 
amounts  of  hysteresis  and  nonlinearity  observed  in  rocks  under  corresponding 
stresses. 

Experimental  results  also  show  that  each  rock  displays  nearly  linear 
behavior  when  the  loading  is  very  small,  around  4.5  x  104  Pa,  resulting  in 
strains  between  2  x  10"®  and  6  x  10"®,  depending  on  rock  type  (Figs.  5.5  through 
5.81.  At  higher  load  levels  (Figs.  5.9  through  5.12)  there  is  evidence  of 
strong  nonlinearity  and  inelasticity  in  the  weakest  frame  rock,  Berea  sandstone, 
while  the  nonlinear  and  inelastic  effects  are  more  subtle  in  the  stronger  frame 
rocks  especially  Indiana  limestone.  In  all  rock  specimens  except  Indiana  lime¬ 
stone  the  effective  modulus  of  the  rock  is  significantly  larger  in  compression 
than  in  tension.  In  one  anomalous  case,  Indiana  limestone  has  a  higher  modulus 
in  tension  when  the  load  is  sufficiently  small  (Fig.  5.13). 


Fig.  5.5  Stress-strai n  response  for  Boise  sandstone  with  a  peak  load¬ 
ing  stress  of  4.5  x  It)*4  Pa  showing  nearly  linear  behavior. 
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srmfll.  Surface  strain 

Fig.  5.6  Stress-strain  response  for  Berea  sandstone  with  a  peak 

loading  stress  of  4.5  *  10H  Pa  showing  nearly  linear  behavior. 


S A M P L t.  SUKr  rtCE  STkRIN 

Fig.  5.7  Stress-strain  response  for  Westerly  granite  with  a  peak 

loading  stress  of  4,5  *  1 0 14  Pa  showing  nearly  linear  behavior. 
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Fig.  5.8  Stress-strain  response  for  Indiana  limestone  with  a  peak 

loading  stress  of  4.5  *  104  Pa  showing  nearly  linear  behavior. 


EPVP’LL'  SL'RFPCE  STRRJN 


Fig.  5.9  Stress-strain  response  for  Boise  sandstone  with  a  peak 
loading  stress  of  1.32  *  106  Pa. 
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Fig.  5.12  Stress-strain  response  for  Indiana  limestone  with  a  peak 
loading  stress  of  1.32  *  106  Pa. 


29 

C7621TC/jbs 


Rockwell  International 

Science  Center 


SC5361.10FR 


;a,'D.E  S_'“  ~  "C  l  E’PP]’; 


Fig.  5.13  Stress-strain  response  for  two  rocks  showing 
the  anomalous  occurrence  of  a  slightly  higher 
effective  modulus  in  tension  than  in  compression. 
Peak  loading  stresses  are  relatively  small 
(1.32  x  105  Pa)  in  each  case. 
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McKavanagh  and  Stacy  (1974)  argue  that  it  is  possible  to  distinguish 
between  linear  anelastic  relaxation  and  nonlinear  relaxation,  such  as  inter¬ 
granular  friction,  by  examining  the  shape  of  the  hysteresis  loop  tips.  Cusped 
tips  indicate  a  nonlinear  mechanism,  while  rounded  tips  indicate  linear  anelas¬ 
tic  relaxation.  Expanded  views  of  the  hysteresis  loop  tips  measured  on  Westerly 
granite,  Boise  sandstone,  and  Bere3  sandstone  are  shown  in  Figs.  5.14,  5.15,  and 
5.16,  respectively.  In  each  case  the  tips  are  cusped  and  not  rounded,  which 
argues  against  an  anelastic  relaxation  mechanism. 
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Fig.  5.14 

Magnified  views  of  the  hysteresis 
loop  tips  for  Westerly  granite. 
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Fig.  5.15 

Magnified  views  of  the  hysteresis 
loop  tips  for  Boise  sandstone. 
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Fig.  5.16 

Magnified  views  of  the  hysteresis 
loop  tips  for  Berea  sandstone. 
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In  some  cases  there  is  clear  evidence  for  large  nonrecoverable  changes 
in  the  rock  during  the  first  excursion  to  a  "new"  maximum  load.  It  is  reason¬ 
able  to  assume  that  these  changes  are  associated  with  intergranular  cracking. 
This  is  apparent  when  the  hysteresis  loop  is  not  clearly  established  until  at 
least  the  second  cycle  of  loading.  In  the  case  of  Berea  sandstone  subjected  to 
1.25  x  10^  Pa  axial  stress  (Fig.  5.17),  the  result  of  the  first  full  loading 
cycle  was  a  net  shortening  of  the  specimen.  The  nonrecoverable  change  at 
comparable  stress  is  small  in 
rocks  with  a  moderately  strong 
frame,  such  as  Boise  sandstone 
and  Westerly  granite  (Figs. 

5.18  and  5.19) .  Indiana 
limestone,  the  rock  with  the 
strongest  frame,  shows  no 
hysteresis  and  no  nonrecover¬ 
able  changes  after  comparable 
loading  (Fig.  5.20). 
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Fig.  5.17 

Comparison  of  the  first  load- 
ingcycle  for  Berea  sandstone 
at  1.25  x  io6  Pa  with  the  sub¬ 
sequent  7  cycles  showing  non¬ 
recoverable  changes  after  the 
first.  The  previous  maximum 
load  was  4.39  *  105  Pa. 
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EOISE  SRNDST ONE  B80-34R 
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Stress-strain  response  for  the  first  eight  loading  cycles  for  Boise 
sandstone  at  1.25  *  IQ6  Pa.  A  small  amount  of  nonrecoverable  change 
is  observed  after  the  first  cycle,  but  the  following  seven  are  super 
imposed.  The  previous  maximum  load  was  4.39  *  105  Pa. 
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Stress-strain  response  for  the  first  eight  loading  cycles  for 
Westerly  granite  at  1.25  *  106  Pa.  A  small  amount  of  nonrecoverable 
change  is  observed  after  the  first  cycle,  but  the  following  seven  are 
superimposed.  The  previous  maximum  load  was  4.39  *  105  Pa. 
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INDIflNR  LIMESTONE  333-4 
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Fig.  5.20  Stress-strain  response  for  the  first  eight  loading  cycles  for  Indiana 
limestone  at  1.25  *  106  Pa.  The  material  appears  to  be  nearly  linear 
elastic  and  nonrecoverable  changes  are  not  observed. 

By  reversing  the  direction  of  initial  loading  it  is  possible  to  examine 
the  shape  of  tensional  hysteresis  loops  and  compressive  hysteresis  loops.  In 
Figs.  5.21  through  5.24  we  show  the  effects  of  reversing  the  loading  direction 
from  tension-first  to  compression-first,  and  from  compression-first  to  tension- 
first,  respectively,  for  two  different  rock  types.  It  is  apparent  that  in  Boise 
sandstone  the  opening  of  the  loop  is  due  primarily  to  tensile  stresses,  and  that 
the  behavior  in  compression  is  nearly  nonlinear  elastic.  In  Westerly  granite 
these  effects  are  similar,  but  the  contrast  between  compressional  hysteresis  and 
tensional  hysteresis  is  less  pronounced.  A  hysteresis  loop  exists  in  both 
compression  and  tension,  but  it  is  smaller  in  compression. 

The  effect  of  frequency  on  the  shape  of  the  hysteresis  loop  has  also 
been  examined.  The  results  of  measurements  on  Boise  sandstone  and  Westerly 
granite  at  1  Hz  and  at  0.1  Hz  are  shown  in  Figs.  5.25  and  5.26.  The  curves 
superimpose  very  well,  which  is  consistent  with  a  frictional  relaxation  mecha¬ 
nism,  even  though  this  evidence  alone  does  not  preclude  anelastic  relaxation. 
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Fig.  5.21 


Effect  of  reversing  the  initial 
loading  direction  from  tension- 
first  to  compression-first  for 
Boise  sandstone  with  a  maximum 
loading  stress  of  4.39  x  106  Pa, 
showing  the  lack  of  significant 
hysteresis  in  the  compressional 
hysteresis  loop. 
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Fig.  5.22 

Effect  of  reversing  the 
initial  loading  direction 
from  compression-first  to 
tension-first  for  Boise 
sandstone  with  a  maximum 
loading  stress  of 
4.39  *  106  Pa,  showing  the 
development  of  a  hysteresis 
loop  in  tension. 
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Fig.  5.23 

Effect  of  reversing  the  initial 
loading  direction  from  tension- 
first  to  compression-first  for 
Westerly  granite  with  a  maximum 
loading  stress  of  4.39  x  106  Pa 
showing  a  significant  com- 
pressional  hysteresis  loop. 


Fig.  5.24 

Effect  of  reversing  the 
initial  loading  direction 
from  compression-first  to 
tension-first  for  Westerly 
granite  with  a  maximum  load¬ 
ing  stress  of  4.39  *  106  Pa. 
Figure  shows  the  existence 
of  a  hysteresis  loop  in 
tension  which  is  larger  than 
the  compressional  hysteresis 
loop  illustrated  in  Fig.  5.23. 
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Fig.  5.25 

Comparison  of  the  stress-strain 
response  for  Berea  sandstone  at 
two  different  frequencies.  The 
traces  appear  to  be  identical. 


Fig.  5.26 

Comparison  of  the  stress- 
strain  response  for  Westerly 
granite  at  two  different 
frequencies.  The  traces 
appear  to  be  identical . 
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Finally,  we  have  estimated  the  effective  attenuation  from  the  areas  under 
the  hysteresis  loops.  In  most  cases  the  hysteresis  loops  used  for  these  calcula¬ 
tions  were  obtained  using  alternating  compressive  and  tensile  loading,  as  shown  in 
Figs.  5.1  through  5.20.  The  results  for  all  four  rock  types  are  plotted  as  a 
function  of  maximum  strain  amplitude  in  Fig.  5.27.  In  two  cases  it  was  also 
possible  to  calculate  Q  from  hysteresis  loops  that  represented  only  tensional  or 
only  compressional  loading,  as  in  Figs.  5.21  through  5.24.  The  results,  shown  in 
Table  5.1,  show  clearly  that  the  effective  attenuation  in  compression  is 
significantly  less  than  in  tension,  although  significant  differences  do  exist 
between  3oise  sandstone  and  Westerly  granite. 

ATTENUATION  VS  AMPLITUDE 


Fig.  5.27 

Attenuation  calculated  from  whole 
cycle  hysteresis  loop  areas  as  a 
function  of  maximum  strain  amplitude. 
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Table  5.1 

Estimates  of  Q  from  Hysteresis  Loop  Areas  Measured  in  Tension  Only, 
in  Compression  Only,  and  in  Both  Tension  and  Compression. 


TENSION 

COMPRESSION 

TENSION  &  COMPRESSION 

WESTERLY  GRANITE 
MAXIMUM  STRAIN  -  1 .5  E-4 

186 

23  6 

21  2 

BOISE  SANDSTONE 

MAXIMUM  STRAIN  .  2  E-4 

22.5 

200 

41.8 
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6.0  CONCLUSIONS 

We  have  evaluated  the  techniques  currently  in  use  in  the  laboratory  for 
measuring  the  attenuation  of  seismic  waves.  We  conclude  that  the  forced  swept 
resonance  technique  is  best  suited  for  measuring  the  Q  of  linear  anelastic  ma¬ 
terials,  and  for  defining  the  transition  amplitude  from  linear  anelastic  behavior 
to  nonlinear  behavior.  For  evaluating  losses  in  materials  at  high  amplitudes, 
which  exceed  the  elastic  limit,  the  best  technique  requires  simultaneous 
measurements  of  the  time  functions  of  stress  and  strain  under  conditions  that 
simulate  seismic  loading  as  closely  as  possible. 

Using  the  mechanical  resonance  approach  we  have  studied  the  transition 
region  from  linearity  to  nonlinearity  in  a  test  specimen  of  Westerly  granite  at 
elevated  effective  pressures.  Nonlinear  effects  in  shear  are  observed  when  the 
maximum  shear  strain  exceeds  approximately  10‘6,  increasing  slightly  with  in¬ 
creasing  effective  pressure.  Nonlinear  effects  in  flexure  are  also  observed  when 
the  strain  exceeds  10-6,  also  increasing  very  slightly  with  increasing  effective 
pressure.  These  transition  amplitudes  probably  represent  a  lower  limit  on  the 
amplitude  of  transition  from  linear  to  nonlinear  behavior  for  the  compressional 
pulse  propagating  through  the  near-field  of  an  explosion. 

In  this  study  we  have  also  examined  experimentally  the  details  of  stress- 
strain  hysteresis  loops  when  a  number  of  different  rock  specimens  are  subjected  to 
alternating  compressive  and  tensile  stresses.  In  general  at  nonlinear  amplitudes 
the  stiffness  of  the  rock  is  greater  in  compression  than  in  tension.  Furthermore, 
preliminary  results  indicate  that  most  of  the  energy  loss  during  a  full  cycle  of 
loading  occurs  as  a  result  of  strain  in  extension:  the  hysteresis  loop  in  compres¬ 
sion  is  smaller  than  the  hysteresis  loop  in  tension.  The  shape  of  the  hysteresis 
loops  also  appear  to  be  independent  of  frequency.  These  observations  indicate  a 
loss  mechanism  associated  with  intergranul ar  friction.  Intergranul ar  sliding 
appears  to  be  restricted  by  the  impingement  of  opposing  crack  faces  in  compres¬ 
sion,  since  the  rock  is  stiffer  in  compression  than  in  tension  and  a  large 
hysteresis  loop  develops  only  when  the  rock  is  subjected  to  tensile  stresses. 
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The  results  of  this  study  indicate  that  the  mechanical  behavior  of  rocks 
can  be  significantly  different  in  compression  than  in  tension  and  that  the  onset 
of  nonlinear  effects  with  increasing  strain  may  not  be  the  same  for  tensile  loads 
as  for  compressive  loads.  All  available  evidence  indicates  that  the  primary 
relaxation  mechanism  at  nonlinear  amplitudes  between  10“®  strain  and  10'^  strain 
involves  i ntergranul ar  friction.  More  experimental  work  in  this  area  will  shed 
light  on  the  issue  of  linearity  vs  nonlinearity  at  intermediate  strains,  and  also 
will  provide  realistic  detailed  information  about  rock  rheology  for  the  numerical 
modeling  of  near-field  seismic  pulse  propagation.  For  this  work  to  be  most 
meaningful  true  uniaxial  strain  is  necessary,  and  specimens  must  be  exposed  to 
elevated  confining  pressures. 
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